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Abstract 

We propose a mechanism of thermalization of nucleons in relativistic nucleus-nucleus collisions. 
Our model belongs, to a certain degree, to the transport ones; we consider the evolution of the 
system, but we parametrize this development by the number of collisions of every particle in 
the system rather than by the time variable. We based on the assumption that the nucleon mo- 
mentum transfer after several nucleon-nucleon (-hadron) collisions becomes a random quantity 
driven by a proper distribution. 
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The Model. Our model is aimed at the description of the nucleon spectra {dN/dy, 
dN/m±dm±dy) for such colhsion energies, when the number of created nucleon-antinucleon 
pairs is much less than the number of net nucleons, i.e. it can be applied at AGS and 
low SPS energies. 

We now make three key assumptions about the nucleon system. 

1. We separate all nucleons after freeze-out into two groups in accordance with their 
origination: a) the first group consists of net nucleons that went through hadron reactions; 
b) the second group includes nucleons which were created in the collective processes, for 
instance, during the hadronization of a QGP. Then, the nucleon momentum spectrum 
can be represented as a sum of two different contributions: 

dN fdN\ {'^^^ 



d^p Vc^^^y hadron ' K^^P J QGP ' 

In turn, the total number of registered nucleons equals iVtotai = -^hadron + Nqqp. If 
this separation can be done, we can define the "nucleon power" of the created QGP as 
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^QGP ~ -^QGP /-^totai ■ In the present work, we mainly deal with the nucleons from the 
first group. 

2. The collision number for every nucleon (hadron) is finite because the lifetime of the 
fireball is limited. To determine the maximal number of collisions, Mmax, in a particular 
experiment, we use the results of UrQMD calculations |l|2j . 

3. Because the colliding nuclei are spatially restricted many-nucleons systems, the 
different nucleons experience different collision numbers: it is intuitively clear that the 
collision histories of the inner and surface nucleons will be different. That is why, we 
partition all amount of nucleons of the first group (nucleons which take part in hadron 
reactions only) into different ensembles in accordance with a number of collisions before 
freeze-out. Then the nucleons from every ensemble give their own contribution to the 
total nucleon spectrum. If we denote the number of particles in a particular ensemble 
where the particles experienced M effective collisions by C(M), then, in correspondence 
to ll]), we can write the total nucleon spectrum as 

A/max 

— = ^ C{M) Dm{p) + Cthcrm Ahcrm(p) , (2) 

^ M=l 

where Dm{p) is the spectrum (normalized to unity) of the particles in the A/-th ensemble. 
The last term on the r.h.s. of ([2]) corresponds to the possible contribution from the 
totally thermalized source which we associate with a QGP. Here, Cthcrm is the number 
of nucleons which are created during the hadronization of the QGP, and -Dthcrm(p) is the 
thermal distribution. 

Consider successive variations of the momentum of the n-th nucleon from nucleus A 
which moves with momentum fep along the collision axis from left to right toward nucleus 
B. Every m-th collision induces the momentum transfer ql™^ for the n-th nucleon. So 
that, after M collisions, the nucleon acquires the momentum fc„ = feg + Qn^ where 
= X]m=i 9""'' is total momentum transfer finally obtained by the n-th nucleon. 

We assume for the moment that the elastic scattering gives the main contribution to 
the two-nucleon collision amplitude. The initial momentum of every nucleon in nucleus 
A is ka = — (0,0, /cqz), while the initial momentum of every nucleon in nucleus B 
is fcb = — feo = (OjOj ~^Oz)i in the c.m.s. of colliding nuclei. The energy and momentum 
are conserved in every separate collision of two particles, Lo{ka) + u){kij) — uj{Pa) + 
uj{pi,), ka + kb — Pa + Pbi where p^ and are the momenta of the particles after 
collision. We assume that the particles are on the mass shell, so that uj{k) — \/ irfi + k^ 
(the system of units ?i = c = 1 is adopted). Determining the six unknown quantities, p^ 
and pj, from four equations is straightforward but two quantities (e.g., {Pa)x and {Pi,)x) 
remain uncertain and can be considered as such which accept random values driven 
by the scattering probability. After the third collision, every component of the particle 
momentum becomes completely random. If the initial momentum is fixed, this means 
the full randomization of the momentum transfer after three successive collisions. So, if 
we follow the elastic scattering of a nucleon from the first collision to the last one, we 
would see the full randomization of the momentum transfer after every three successive 
scattering acts. As for the inelastic collisions, the nucleon momentum transfer undergoes 
the even faster randomization [3j . 

First, we would like to determine the density distribution function f2N in the mo- 
mentum space, which describes 2A'^ nucleons after M collisions per particle. The whole 
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consideration is carried out in the c.m.s. of two identical^coUiding nuclei. Let us write 
down a density distribution function in the form f2N = Cf2N, where C is the normahza- 
tion constant. The unnormaUzed distribution function f2N can be defined in a two-fold 
way: first, we follow all collisions of a particular nucleon by the integration with respect 
to all nucleon random momentum transfer, and, second, we fix the total energy of the 
2iV-nucleon system after freeze-out in a microcanonical-like way: Etot = Sn=i'^(^n)- 
The integration measure of the momentum transfer for the n-th particle in a series of M 
collisions is determined as 



M ,3 (m) n 13 



(3) 



where the distribution of the momentum transfer is characterized by the presence of the 
form- factor Jm(<7). For the sake of simplicity, we assume the independence of Jmio) on 
the collision number m, i.e. Jm{q) ^ J{q)- Hence, we adopt an approximation where just 
one form-factor J{q) characterizes a distribution of the momentum-transfer in a series 
of collisions which are experienced by a nucleon during its traveling through the fireball. 
Then the unnormalized 27V-particle distribution function reads 



f2N{Etot', ki, 

N 



2N 



n 



(27r)-^(5^ k 




M 

E 



,(m) 



(4) 



where V is the volume of the system in the coordinate space. Making the Laplace trans- 
formation, we determine the unnormalized distribution function 



F2^(/3;fci 



27V) 



dEtote~^^'°'f2N{Etouki,...,k 



2N) 



(5) 



The distribution function in the canonical ensemble reads ¥2N{f3) = ¥2N{f3)/ Z2n(P), 
where Z2n{P) is the partition function. Skipping over details, we write the result 

N 2N 



=N+l 



where 

fa{b){k) 



-liu)(k) 



(fik 

(2^ 



-I3u}(k) 



iMikTko) (7) 



are the single-particle distribution functions attributed to nucleus "A" for subindex a or 
to nucleus "B" for subindex b, respectively. Note that Za{l3) — zii{/3) = z{(3) for identical 
nuclei. We define the "multi-scattering form-factor" 



Tm{Q) ^yj d're-^Q- [J{r)]'' with J(r) = J A^J{q)e^i-r, 



(8) 



For a large enough number M of the effective collisions and for the form-factor which 
possesses the spherical symmetry, J{q) = J{\q\), we can calculate the first integral in ^ 
within the saddle-point method 
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Im{Q) 



3/2 



where {q^ 



(9) 



^M(g2) J - ' ' ' J (27r)3 

In the limit case M — > oo, the dependence on the initial momentum fcg is washed out, 
and both single-particle distributions fa(b) C^) take the same "thermal" limit: fa(b) (^) ~^ 

/thcrm(fe) = e-'^'^('=)/zthcrm(/3), whcrC Ztherm(/3) = / d^fc e-'^^C") . 

Using a "two-source" single-particle distribution function /(fca, fcfe) = fa{ka) fb{kb) to 
average the quantity W{ka, kb) = \ \5^{p — ka) + 5^{p — kb)\ , we obtain 



Dm{p) 



1 dN\ 



e-^^^P^ [lMip-ko) + lMip + ko)]. 



(10) 



It is evident that the spectrum has two items which can be attributed to the first and 
second colliding nuclei, respectively, and hence it can be named a "two-source single- 
particle spectrum" . 

Nucleon Rapidity Distribution and Transverse Spectrum. To obtain the trans- 
verse mass and rapidity distributions, we pass to new variables: mj^ = {m^ + p^Y^"^ , 
p\ ^ Px ~^ Pyi tanhy = p^/tOp, then d^p = d(j)ujp m± d'm± dy. In accordance with 



d^N 



m±dm±dy 



= 27r m± coshy 



J2 C{M)DMim^,y)+ Cth crm Ah crm 



(11) 



We can define the distribution functions in the rapidity space as 

^M{y)= TTCOshy dm±ni\[fa{in±,y) + fb(in±,y)] with / $m(2/) = 1 ■ (12) 
Then, the rapidity distribution looks like 



dN 
dy 



M„- 



C{M) ^M{y) + Cthcrm «'thcrm(y) 



(13) 



M=l 



Toy Model. To present the explicit results of our approach, we consider a toy model: 
the form-factor J(g) is chosen as a homogeneous distribution in the sphere of finite radius 



(2^)^ 
V„ 
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with 



d3 
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(27r)3 



1. 
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In the proposed model, the maximum number of collisions, Mmax, is assumed to be 
finite and determined by the nucleus number initial energy, and centrality. With the 
help of the UrQMD transport model [1,2J, it was found that, under the AGS conditions [4] 
with a centrality of (0-3)%, M^ax = 13 and Qmax — 0.8 GeV/c (for the toy model | q^^^ = 
(q^) = (9^) UrQMD)- Utilizing the thermal distribution, we extract a slope parameter from 
experimental data on the proton m^-spectra [4], T = 280 MeV. Note, the proton data are 
of interest, first of all, because we know an exact value of the initial nucleon momentum. 
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AGS(E802): Au+Au, p =11.6AGeV/c 



AGS(E802): Au+Au, p,_j=11.6A GeV/c 
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Fig. 1. The result of a fit (thick sohd curve) to 
experimental data [4] on the rapidity distribution. 
Broken curves marked by numbers M represent 
the partial contributions from every collision en- 
semble, C{M) ^M(y) ■ 




Fig. 2. Solid curves represent mj^-spectra ob- 
tained in accordance with JTTJ. We use values of 
the coefficients C{M) which were obtained as a 
result of the fit to the dN/ dy data. Experimental 
points are from W . 



The results of the fit to experimental data |4j on the rapidity distribution and mj_- 
spectra of the net protons are depicted in Figs. [H [51 The fit was done with making use of 
expansion (|f 3p and resulted in a good description of the experimental data. It turns out 
that, in the case of a small number of experimental points, the set of functions <&m(2/) is 
overcomplete. To choose a unique configuration of the variable parameters C(M), we use 
the maximum entropy method [5]. The set of coefficients C{M) (see Table I) are nothing 
more as the absolute number of protons in every collision ensemble. This description 
was carried out with account for the contribution of the thermal source, which cannot 
appear due to the nucleon-nucleon (-hadron) rescattering (see the partial contributions 
in Fig. 1). We assume that this source is a thermalized multiparton system (QGP) which 
emits totally thermalized nucleons through the hadronization process. The knowledge 
of the number of protons, Cthorm, which come from the QGP, gives us a possibility to 
evaluate the "nucleon power" of the QGP, -Pqgp, created in a particular experiment on 
the nucleus-nucleus collision. We find that, under the AGS conditions [1] (a centrality of 
0-3%) 

' -^QGP ~ So, in the framework of the proposed criterion, it could be claimed 
that the QGP (as a nucleon source) was created not only at SPS energies p|, but it was 
also created in the central collisions at AGS energies. 

Meanwhile, the partial expansion, dN/dy = X)m=T C'(M) <pM{y) (see (dS])), where we 
do not use the thermal contribution, makes a good description (the same x^) of the ex- 
perimental data on both the rapidity distribution and m_L-spectra. So, we cannot resolve 
unambiguously the presence of the thermal source. In fact, to overcome the problem, we 
need a more detailed experimental information for the central rapidity region. 

All this encourages us to apply the model to other experiments and problems. 
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